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Abstract 

A model magnetosphere for Mercury is presented using an upstream 
image-dipole and nightside 2-dimensional tail current sheet method. The 
tail field is represented by a new analytical formulation. The magnetic 
field data from the Mercury I encounter by Mariner 10 in March 1974 are 
used to determine quantitative parameters of the model magnetosphere, 
using the method of least squares. The magnetopause crossing points 
directly observed are used to determine the size of the magnetosphere, and 
the solarwind conditions are used to determine the magnetospheric field 
at the stagnation point. The model produces a magnetosphere- like region 
with planetary field lines that are confined in nearly circular cross- 
sections transverse to the Sun-planet line : Results ire used to show 

geometry, field line configuration, and contours of constant field 
intensity inside the magnetosphere. 



Introduction 


This paper presents a quantitative model of the magnet ospheric 
field of Mercury based on the magnetic field data observed during the 
first encounter of the Mariner 10 spacecraft with Mercury on 29 March 
1974 (hereafter referred to as Mercury I) . 

The preliminary results from the analysis of the Mercury I data 
(Ness et al . , 1974) unexpectedly revealed that the planet Mercury 
appeared to have a modest but significant intrinsic magnetic field. The 
planetary field was compressed and confined by the dynamic pressure of 
the solar wind plasma to form a magnetosphere- like region. The inter- 
action of the solar wind with Mercury appears Earth-like and includes 
the formation of a detached bow shock upstream of the planet and magnetic 
tail downstream. Simultaneous measurements of the low energy electron 
flux from Mercury I by Ogilvie et al (1974) provided strong correlative 
evidence for the existence of a planetary magnetosphere and a bow shock. 
In addition, intense bursts of high energy charged particles were 
reported on the same spacecraft (Simpson et al. 1974) as occurring in 
a region of space corresponding to the magnetosphere and magnetosheath. 
The. preliminary analysis of Mercury I data yielding an estimate of the 
intrinsic planetary field was based on a simple offset tilted dipole 
model which did not consider external sources of the field. 

Important features of the magnetospheric field of Mercury were 
established in a recent paper by Ness et al. (1975a) in which the field 
was represented by a set of spherical harmonic functions including 
external sources. The planet Mercury occupies a very large fraction 
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of the volume of the magnetosphere since the sub solar point of the 

magnetopause is only ~0„6 R above the surface of the planet. Thus a 

substantial part of the observed magnetosphere field, S o ^ g , is due to 
external sources and net only the intrinsic field. The final analysis showed 

that the centered planetary dipole is nearly normal (83°) to the equatorial 
plane with a polarity sense similar to the Earth's, Significant distortion 
of the observed magnetospheric field from a dipole field was interpreted 
as due to the compression by the solar wind and the formation of a 
magnetic tail. The results suggested that the tail current sheet must 
be located close to the planet on its dark side. 

A second encounter (Mercury II) by Mariner 10 occurred on 21 September 
1974. The targeting strategy of the Mercury II mission was to provide 
optimum imaging coverage of the south polar regions and thus the space- 
craft did not approach sufficiently close to the planet to directly 
observe the magnetospheric field or bow shock of the planet. A third 
encounter (Mercury III) occurred on 16 March 1975 and additional 
observations of the magnetospheric field of Mercury were performed 
which confirmed the existence of an intrinsic field of the planet 
(Ness et al . , 1975b). 

Model Magnetosphere 

This model study is an attempt to determine some quantitative 

parameters for the magnetospheric field of Mercury based on Mercury I 

data. The observed magnetospheric field data available from Mercury I 

are restricted to a single pass during a time period of 17 minutes. 

Let F , F , F denote the three components of the observed vector 
x y z 

magnetic field in a Mercury - centered solar Ecliptic coordinate system (ME), 
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where (x,y,z) denote the location of the spacecraft using the planetary 
radius (R- = 2439 km) as a unit of length. 

r 

Because the available data set is very limited, the model magneto- 
sphere under consideration is a very simple one. As shown in Figure 1, 
it is justifiably assumed that the intrinsic planetary dipole moment 

is located at the center of the planet, and is normal to the equatorial 
plane and the solar wind. Ness et al « (1975a) have determined 
the orientation of the moment to be 0 = -80° and 0 = 285° in the ME 
coordinate system so that the "aspect” angle of solar wind flow relative 
to the dipole axis is 87. 4°, very close to the 90° assumed. 

The tail current sheet is assumed to coincide with the equatorial 
plane on the dark side of Mercury with its edge located at x = x^. We 
assume that the magnetospheric field is divergent- free and curl-free 
everywhere except at the tail current sheet and that the field lines 
are tangential to a magnetosphere- like boundary, i.e. the magnetopause, 
so that the surface is modeled as a tangential discontinuity. 

During the inbound crossing of the magnetopause of Mercury I 
encounter, the spacecraft was located at 


and 


x = -1.66 


(y 2 + z 2 ) 1/2 = 2.30 


The outbound crossing occurred at 
x = -0.78 

and 

(y 2 + z 2 ) 1/2 = 2.30 


Thus the observed magnetosphere during Mercury I encounter had a radius 
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of approximately 2.3 between x = -0,8 to -1.7. The stand-off distance 
to the subsolar point of the magnetopause was not directly observed. 
However, from the solar wind flow velocity and density measured by the 
plasma science experiment (Ogilvie et al 0 , 1974), the magnitude of the 
equivalent magnetic field at the stagnation point to balance the solar 
wind momentum flux is computed as « 160 gammas. 

These conditions will be used in determining quantitative parameters 
for a model magnetosphere of Mercury. Two different models shall be 
developed: (1) assumes the magnitude of the planetary dipole moment to be 

22 3 

unknown, and (2) constrains the moment to the 5.1x10 Gauss-cm value 
determined by Ness et al. (1975a). The latter case corresponds to an 
equatorial field strength of 350/ for the internal dipole. 


Tail Field 


The tail field is represented by the simple analytical form 


B --rising B, 

X t z t 


« - r t ' 1/2 cosf B t 


in the region: r > 0 and -it < to < tt 

HerS [, ,2 2 l 1/2 

r = (x - x t ) 2 + z- 

and ^ ^ 

co = tan 1 ( — ) 

x - x t J 

where r is the distance from the edge of the two-dimensional current 
t 

sheet. The strength of the tail field is characterized by a constant 
B t which represents the magnitude of the tail field at r^= 1. The 
component s 0 since the current sheet extends to infinity in both 
the + and - y directions. 

Equation (1) can be rewritten in the following form as a function 
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of a complex variable 


B + iB = (r 
z x t t 


( 3 ) 


where is a real constant. This function is single-valued and analytic 
at all points in the region > 0 and -tt < ou < rr. Satisfaction of the 

Cauchy-Riemann conditions shows that the tail field given by (1) is 
divergent- free and curl- free everywhere except along a branch cut 
defined by z - 0 and x ^ x which corresponds to an infinitely thin tail 
current sheet. 

In the absence of other fields , the tail field represented by (X) 
immediately on each side of the current sheet is parallel to the sheet 
but in opposite directions on opposite sides. The potential lines are 
represented by 

1/2 . ti> 

r t sm — = constant 


and the magnetic field lines by 

1/2 a) . _ 

r cos ^ = constant 

Figure 2 shows the geometry of magnetic field lines of the tail current 
sheet. 

Magnetospheric Field 

The image-dipole method (Hones, 1963; Taylor and Hones, 1965) 
is used to represent the contribution to the magnetospheric field of 
Mercury from the electrical currents associated with the deflected solar 

wind flow, which forms the magnetopause. Let 

r* o-> 

M = g- e 

1 z 

denote the moment of the centered planetary dipole. An image dipole 
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with a moment of 

3. = r3 

i 

is placed at (x^, 0, 0), in the sunward direction. The image dipole 
represents the effects of confining the planetary field within a magneto- 
pause. Adding the tail field (1) to the field due to M and M^ } we write 
the net field vector for the model magnetosphere of Mercury as 

B =-g° (c e + c e + c e)+B^(c e w + c ej (6) 

model I i I x 21 V 31 2 1 2 x 32 2 

where 



Here r denotes the distance from the center of Mercury, and r_^ the 
distance from the image dipole. The model field (6) is tangential to a 
boundary surface which corresponds to the magnetopause, and it is 
divergent-free and curl-free everywhere within the magnetosphere except 
at the tail current sheet. 

Method of Least Squares 

For a given set of N data points , the vector field residual, or 
mean squared deviation between theory and observation can be written as 
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° 2 - ~ E(Wi - s obs ) 2 

1 n=l 


(7) 


N r 


- V (c g° + c B - F ) + (c g°-F) + (c g° + c B - F ) 

N-l ^ L H I 12 t x' 21 1 y 31 1 32 t z' 


2 . 


is controlled by five parameters: g^, B fc , F, x_^ and x fc0 The minimum 


2 . 


of c is the best fit to the data and is determined by setting — — = 0 

Sg 1 


and = 0 


Thus, a system of two normal equations results: 


b g + b B = b 
i 1 1 i2 t 13 


(i = 1,2) 


( 8 ) 


where 


b = £ (c 2 + c 2 + c 2 ) 


1 1 


n=l 


1 ). 21 


31 


b = b = V (c c + c c ) 
12 21 ^ 1112 31 32 


N 

b = V ( C F + c F + c F ) 
13 11 X 21 Y 31 Z 


b = 52 (c 2 + c 2 ) 

22 n== ] 1 2 32 


N 

b = S (c F + c F ) 

2 3 n =■ | 1 2 X 32 Z 

Making use of (7) and (8), for any given set of N data, B t> g°, (if 
2 

necessary), and a as functions of the three model parameters T, x., 

and x,_ can be calculated, 
t 

The best fit model magnetosphere is thus governed by these three 
parameters: F, the ratio of the moment of the image dipole to that of 



- 9 - 

Mercury's dipole, the distance between the two dipoles and x the 
position for the edge of the tail current sheet „ It is necessary to 
adjust these three parameters until a best-fit model magnetosphere is 
obtained. The geometry of the magnetopause is mainly controlled by 

r and x^. for any given pair of T and x^, the value of x^ is adjusted 

2 

until a minimum value of a is obtained, 

Magnetosphere with Unknown Planetary Dipole Moment (Model 1) 

The Mercury I magnetic field data during the intervals from 2038 
to 2047 UT and from 2049 to 2052 UT (Ness et al } 1975a) have been used 
to construct a representative magnetospheric field of Mercury. When the 
spacecraft was inside the magnetosphere of Mercury, abrupt increases in 
the flux of charged particles have been observed (Simpson et al„ 1974). 

Inside the magnetosphere this observed feature at ~g048 UT was designated 
as the B event and that at ^053 UT as the C event. (The D event occurred 
outside the magnetosphere and the A event was quite small in comparison). 
Magnetic field data during the B and C events are not used in this study. 

To meet the condition that the net magnetospheric field is approximately 
160 gammas at the stagnation point, a series of possible solutions was found 
which corresponded to a varying value of P. The size of the 
magnetosphere increases as the dipole ratio f decreases. The geometry 
of the magnetosphere that best-fits the magnetopause crossing points 
directly observed during Mercury i for the case of g° variable has a 
value of cr = 16 „ 5 gammas. At the radial distance to the sub-solar 
point of the magnetopause, x sta g == 1-33, the dipole field is "bmplif ied" 
by a factor of 2.35. A comparison of the model magnetospheric field 
with the observed field is shown in Figure 3. The value of g° gives 


a 
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dipole moment of the intrinsic field = 2.3x10 Gauss -cm , which is 

approximately one half that obtained from the earlier studies. The 
other pertinent parameters are x.= 6.44, V = 93.3, x = -1.5 and 
B = -35.1 gammas. 

Figure 4a shows the geometry of the model magnetopause and 
the field lines within the magnetosphere. The cross-sections of the 
magnetopause on planes normal to the Sun-Mercury line are nearly circular. 
Their radius near x = - 1 is approximately 2.3 which is in good agreement 
with the observed crossing points of the magnetopause. 

Figure 5a shows the iso- intensity field contours of the model 
magnetosphere in the noon-midnigh t meridian plane. The presence of the tail 
current sheet introduces a neutral line just forward of the sheet on 
the equatorial plane. The singularity line along the edge of the sheet 
is a mathematical construct due to the assumed infinitely thin tail current 
sheet in this model. 

Magnetosphere with Known Planetary Dipole Moment (Model 2) 

In case the planetary dipole moment is known, the method of least 
squares leads to a single normal equation which determines the 
characteristic magnitude of the tail field 

B t = (b S3 +b 21 g l )b 22 (9) 

Assuming that g° =-350 gammas for Mercury, the same set of Mercury I 
magnetic field data is used to study another model of the magnetosphere 
of Mercury. To meet the condition that the magneto spheric field is 
approximately 160 gammas at the stagnation point, a series of possible 
solutions corresponding to a varying value of T is again obtained. 

Model magnetospheres determined by this series of solutions have 
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cr-values always greater than 37 gammas. Figure 4b shows the geometry 
of the magnetopause and field lines in the noon-midnight meridian plane 
represented by a model with fixed g° =-350 gammas . Figure 5b presents 
the corresponding iso- intensity field contours in the noon-midnight 
meridian plane. In this case, the stagnation, point distance is found 
to be 1.71 with a tail field B = -91.67 and = -1.95 while the image 
dipole is located at = 10.4 with a strength F = 264 which combine to 
yield a = 46/. The fit to observations is shown in Figure 3. 

Discussion 

Quantitative models of the magnetosphere of Mercury have been 
developed which are useful in studying the interaction of the solar 
wind with Mercury and the magnetic field geometry within the 
magnetosphere. When the momentum flux of the solar wind changes, the 
stand-off distance and the magnetospheric field at the stagnation point 
varies. The present models predict that the stand-off distance will vary 
between 1.2 and 1.7 under varying solar wind conditions and that 
theoretically it is very unlikely that the solar wind can directly 
impact the surface of Mercury except via access at neutral points 
in the polar regions. This problem has also recently been studied 
empirically by Siscoe and Christopher (1975), who reached the same 
conclusion using solar wind measurements at 1 AU, scaled to the orbit 
of Mercury, i.e. between 0.31 - 0.47 AU. 

The present model wilL be best suited for study of the magneto- 
spheric field near the noon-midnight meridian plane and will be less 
valuable when approaching the downstream flanks. This is due to the 
use of a 2 dimensional tail current sheet, rather than a restricted 
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3 dimensional tail current sheet. Constraining the dipole term leads 
to a substantially inferior fit to the observed data at Mercury I. 
Although the direction of the field is moderately close, the magnitude 
is consistently larger than actually observed. 

The observations at Mercury III will help in ascertaining which of 

o 

the two models, and hence g-^ values, are closer to the actual case. It 
may also be expected that relaxing the condition that the observed 
magnetopause crossings agree well with the theoretical model will 
yield a more realistic field configuration in the region near the 
noon-midnight meridian plane, |y| <1.5. This will be especially true 
if the stagnation point distance is predetermined by a scaling from an 
analogy with the terrestrial magnetosphere. 

In either of the models, there is little if any indication of the 
development of a field configuration which could confine charged 
particles in a "trapped” or radiation belt structure. This is because 
of the large day-night asymmetrices of the field associated with its 
^ i §ti.i £ xcan t compression by the solar wind on the dayside and extension 
on the nightside to form a magnetic tail. The magnetospheric field is 
so distorted that a careful study of an accurate model shall be required 
to further quantify these remarks. There are large spatial gradients 
of the magnetic field intensity so that slight changes in the solar wind 
momentum flux and direction, lead to significant changes in the 

comparisons with actual observations, A more accurate model would then 
begin by using a solar magnetospheric coordinate system to account for 
the 10° obliquity of the dipole axis to the ecliptic and a 3 dimensional 
confined tail current sheet. 
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Figure 4 


A model magnetosphere of Mercury includes a tail current 
sheet coincide with the equatorial plane on the darkside of 
the planet* The tail current sheet has a straight 
edge located at x = x . 

The field line topology of the 2 dimensional tail field. 

The field intensity F, latitude 9 and longitude 0 of the 
model magnetospheres compared with observed data along the 
trajectory of Mercury I. 

Field lines and magnetopause shapes in the noon -midnight 
meridian plane of the model magnetospheres of Mercury. 

Field iso-intensity contours in the model magnetospheres 
in the noon-midnight meridian plane. 
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MODEL MAGNETOSPHERE OF MERCURY 
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FIELD LINES OF MODEL TAIL CURRENT 
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